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In this paper we study the cosmological dynamics of Randall - Sundrum braneworld type scenarios 
in which the five - dimensional Weyl tensor has a non - vanishing projection onto the three-brane 
where matter fields are confined. Using dynamical systems techniques, we study how the state 
space of Friedmann-Lemaitre-Robertson-Walker (FLRW) and Bianchi type I scalar field models 
with an exponential potential is affected by the bulk Weyl tensor, focusing on the differences that 
appear with respect to standard general relativity and also Randall - Sundrum cosmological scenarios 
without the Weyl tensor contribution. 

PACS numbers: 


I. INTRODUCTION 

The notion that we live on a three - dimensional brane 
embedded in a higher - dimensional spacetime has at¬ 
tracted a considerable amount of interest over the last 
few years. These ideas have had a long history (for refer¬ 
ences see |jj) but its recent revival is due to recent work 
by Randall and Sundrum j|]. 

In such braneworld scenarios the ordinary matter is 
confined to the brane while gravity can propagate in the 
whole spacetime (known as the bulk). The effective four - 
dimensional gravity on the brane is modified by extra 
terms ]|], Q] in the gravitational equations, one quadratic 
in the energy - momentum tensor and the other represent¬ 
ing the electric part of the five - dimensional Weyl tensor. 

The dynamics of a braneworlds filled with a perfect 
fluid have been intensively investigated during last three 
years. In particular Campus and Sopuerta i § have 
recently studied the complete dynamics of Friedmann- 
Lemaitre-Robertson - Walker (FLRW) and the Bianchi 
type I and V cosmological models with a barotropic equa¬ 
tion of state using an approach first introduced by by 
Goliath and Ellis JtJ (see also || for a detailed discus¬ 
sion of the application of dynamical systems techniques 
to cosmology). 

Their analysis led to the discovery of new critical points 
corresponding to the Binetruy - Deffayet - Langlois (BDL) 
models Q|, representing the evolution at very high en¬ 
ergies, where effects due to the extra dimension become 
dominant. These solutions appear to be a generic feature 
of the state space of more general cosmological models. 
They also showed that the state space contains new bi¬ 
furcations, demonstrating how the dynamical character 
of some of the critical points changes relative to the gen¬ 
eral-relativistic case. Finally, they showed that for mod¬ 
els satisfying all the ordinary energy conditions, causality 
requirements and for gamma > 1 , the anisotropy is neg¬ 
ligible near the initial singularity, a result first demon¬ 
strated by Maartens et. al. ||. 

Scalar field dynamics on the brane is considerably more 
interesting due to the fact that the equation of state 
parameter 7 is now dynamical. In a recent paper fic|l 


we considered the dynamics of inflationary models with 
an exponential potential V(<f> ) = exp(b(f>), an important 
class of inflationary models first considered by Burd and 
Barrow |ll| and Haliwell [|l2| in the context of standard 
general relativity (GR). These models have a richer dy¬ 
namical structure in the brane - world scenario and have 
the nice feature that inflation can take place with poten¬ 
tials ordinarily too steep to sustain inflation which re¬ 
sult from the high-energy corrections to the Friedmann 
equation |^3|. Indeed for b < 0 we found past attrac¬ 
tors representing high energy steep inflationary models 
in which inflation ends naturally as the energy drops be¬ 
low the brane tension and the condition for inflation no 
longer holds. 

In this paper we complete this analysis (hereafter ref¬ 
ereed to as Paper I), studying the effects of the five- 
dimensional Weyl tensor on FLRW and Bianchi type 
I cosmological models. In the case of FLRW models 
this study will be completely general whereas in the 
Bianchi type I case we will neglect (or constrain) the 
Weyl tensor components for which the theory does not 
provide evolution equations. The paper is organized as 
follows. In section [n] we briefly outline the geomet¬ 
ric formulation of brane - world scenarios and give the 
main dynamical equations for FLRW and Bianchi type 
I cosmological models where the matter is described by 
a dynamical scalar field <fi with an exponential poten¬ 
tial V((f>) = exp (b<j>). The dynamics of the FLRW and 
Bianchi type I cosmological models are presented in sec¬ 
tion III and IV respectively. Finally we end with a dis¬ 
cussion of the main results of the analysis in section 
Through this paper we follow the following notation: up¬ 
per - case Latin indices denote coordinates in the bulk 
spacetime ( A , B,... = 0 — 4) whereas lower - case Latin 
indices denote coordinates on the brane (a, 6,... = 0 — 3). 
We use physical units in which c = 1. 


II. PRELIMINARIES 

To begin with we introduce the geometrical framework 
of braneworld cosmological models and the main assump- 
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tions used to study FLRW and Bianchi type I cosmolog¬ 
ical models. 


A. Basic equations of the brane-world 

In Randall - Sundrum brane - world type scenarios mat¬ 
ter fields are confined in a three - brane embedded in a 
five - dimensional spacetime (bulk). It is assumed that 

(5) 

the metric of this spacetime, g ABl obeys the Einstein 
equations with a negative cosmological constant A( 5 ) 

mm 

Gab = ~^-(5)9ab + k ( 5 )<Kx) [-^ 9ab + Tab] , (1) 

where G AB is the Einstein tensor, K( 5 ) denotes the five- 
dimensional gravitational coupling constant and T AB 
represents the energy - momentum tensor of the matter 
with the Dirac delta function reflecting the fact that mat¬ 
ter is confined to the spacelike hypersurface x A = % = 0 
(the brane) with induced metric g A B and tension A. 

Using the Gauss - Codacci equations, the Israel junc¬ 
tion conditions and the Zi symmetry with respect to the 
brane the effective Einstein equations on the brane are 

G a b = A.gab T tC T a b T ab Sab i ( 2 ) 

where G a b is the Einstein tensor of the induced metric 
g a b- The four - dimensional gravitational constant n and 
the cosmological constant A can be expressed in terms 
of the fundamental constants in the bulk (k( 5 ), A( 5 )) and 
the brane tension A ] 2 (| . 

As mentioned in the introduction, there are two cor¬ 
rections to the general - relativistic equations. Firstly S a b 
represent corrections quadratic in the matter variables 
due to the form of the Gauss - Codacci equations: 

S ab = ±TT ab - \T a c T bc + ±g ab [3 T cd T cd - T 2 ] . (3) 

Secondly S a b, corresponds to the “electric” part of the 

(5) 

five - dimensional Weyl tensor C ABCD with respect to the 
normals, n A ( n A n A = 1 ), to the hypersurface y = 0 , that 
is 

Sab = C A cBD nC n D , ( 4 ) 

representing the non - local effects from the free gravita¬ 
tional field in the bulk. The modified Einstein equations 
([|) together with the conservation of energy - momentum 
equations \7 a T a b = 0 lead to a constraint on S a b and S a b'. 

V a (Sab-rf 5) Sab)= 0 . (5) 

We can decompose S a b into its ineducable parts relative 
to any timelike observers with 4-velocity u a (u a u a = 
- 1 ): 

Sab = — [( U a Ub + ^ h a bjlA + 2 U( a Qb) + 'Pab\ > 


where 

QaU a = 0 , V(ab)=Vab , P°a = 0 , V ab U b = 0 . (7) 

Here IA has the same form as the energy - momentum ten¬ 
sor of a radiation perfect fluid and for this reason is re¬ 
ferred to as the “dark” energy density of the Weyl fluid. 
Q a is a spatial and V a b is a spatial, symmetric and trace - 
free tensor. Q a and V a b are analogous to the usual energy 
flux vector q a and anisotropic stress tensor 7r a b in General 
Relativity. The constraint equation (|^) leads to evolution 
equations for U and Q a , but not for V a b (see Jisjl). 


B. Scalar field dynamics on the brane 


In this paper we consider both FLRW and Bianchi I 
cosmological models on the brane where the matter is 
described by a dynamical scalar field </> with an expo¬ 
nential potential V(cf>) = exp(&</>). In this case the fluid 
4 - velocity can be written as 

u a = __Y^, u a U a = - 1 , ( 8 ) 

q> 


which makes it automatically orthogonal to the hyper¬ 
surfaces of homogeneity (surfaces of constant </>). With 
this choice of 4 - velocity, the energy - momentum tensor 
T (tj, for a scalar field (f> takes the form of a perfect fluid 
(See page 17 in [^6| for details): 

'Sfia — “t“ ph^a , (9) 


with 


and 


P = (10) 

p = \j )2 - V W ■ (ii) 


where <f> is the momentum density of the scalar field and 
V(4>) is its potential energy. If the scalar field is not 
minimally coupled this simple representation is no longer 
valid, but it is still possible to have an imperfect fluid 
form for the energy - momentum tensor fl7| . 

In the case of FLRW models, the effective Einstein 
Equations (§]) lead to the conditions 

Qa = Vab = 0 , 

and this, through the constraint (|B|), further implies 

D a U = 0 O U = U(t ), (12) 


where D a denotes the covariant derivative associated 
with the induced metric on the hypersurfaces of homo¬ 
geneity (h a b = g a b + u a Ub). The situation is somewhat 
different when one considers Bianchi type I models. This 
time we obtain Q a = 0 but we do not get any restriction 
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on V a b ■ Since there is no way of fixing the dynamics of 
this tensor we will study the particular case in which it 
vanishes and the conditions © also apply. 

The key equations describing the dynamics of these 
models are the modified Raychaudlruri equation 


H 


- H 2 



37-2 
6 
2 U 


K 2 p 



37- 1 P 
3 7 -2 A 


Xk 2 ’ 



(13) 


which is an evolution equation for the Hubble parameter 
= V a it a , and the Friedmann equation constraint (a 
first integral of ©: 


H 2 



_F fl+ i a 2 + i A + ^ i(14) 


where 3 R is the scalar curvature of the hypersurfaces or¬ 
thogonal to the fluid flow and 2 cr 2 = a ab a a b is the shear 
scalar ( a a b = h a c hb d S7 ( c w d ) — Hh a b )• We consider only 
the case of a positive four - dimensional cosmological con¬ 
stant, i.e. A > 0. For Bianchi type I models 3 R vanishes 
whereas for FLRW models it is given by 3 R = 6 ka~ 2 {t). 
On the other hand, the shear vanishes for FLRW models 
and for Bianchi type I models the evolution of the shear 
scalar is 


(a 2 )' = —6Hcr 2 . 


In both cases the evolution equation for U is given by |lf 


U = -4 HU , 


and substituting for p and p from 
energy conservation equation 


( 10 ) and ( |Tl| ) into the 


anisotropic Bianchi I models with W / 0 thus extending 
our recent work which only considered FLRW models 
with vanishing non-local energy density U. The key dif¬ 
ference between this work and a similar study carried out 
by Campos and Sopuerta || is that we have a dynamical 
equation of state parameter 7 which evolves according 
to equation (|lj) above. In the case of an exponential 
potential V{<j>) = exp(&(/>), where b < 0 we obtain 

7 = 3fl7(7-2) + V3^(7-2 )b\f^ ■ (19) 

In order to obtain a compact state space we will re - write 
the above dynamical equations in terms of dimensionless 
coordinates that are appropriately expansion normalized. 


III. ANALYSIS OF FLRW MODELS WITH 
EXPONENTIAL POTENTIALS 

In the first paper of this series jl0| we had to distin¬ 
guish between the case 3 R < 0 ancP’R > 0 when intro¬ 
ducing appropriate expansion normalized coordinates for 
the FLRW models. In this paper, we have to consider 
four different subcases when studying the FLRW mod¬ 
els and each one will have to be normalized by different 
quantities in order to obtain a compact state space: (A) 
U > 0 and 3 R < 0; (B) U > 0 and 3 R > 0; (C) U < 0 and 
3 R < 0; (D) U < 0 and 3 R > 0. The total state space 
is composed of these 4 sectors, which are disconnected 
since trajectories cannot leave one sector and enter an¬ 
other. For a detailed description of the state space, in 
particular the invariant sets and the unusual geometry of 
sector D we refer to jdj]. 


p + Q(p + p) = 0, (15) 


A. U > 0 and 3 R < 0 


leads to the 1+3 form of the Klein - Gordon equation 

ij> + ©<(> + V' (</>)= 0 , (16) 

an exact ordinary differential equation for <j) once the 
potential has been specified. It is convenient to relate p 
and p by the index 7 defined by 


As explained by several authors (see e.g. 10 ) , a com¬ 
pact state space is obtained by using the dimensionless 
general - relativistic variables 


Q =*L 

p ~ 3 H 2 ’ 


1 + — 


3 R 

6H 2 ’ 


Ha = 


A 

3H 2 


( 20 ) 


P = (7 - 1 )P & 7 


P + P 
P 



(17) 


This index would be constant in the case of a simple one - 
component fluid, but in general will vary with time in the 
case of a scalar field: 


7 = 0 7 ( 7 - 2 )- 27 X- . (18) 

9 


Notice that this equation is well-defined even for 4> —> 0, 
since Z = i-. 


<P p 

In the next two sections we study in detail the dy¬ 
namics of (a) FLRW models with U ^ 0 and (b) the 


together with the following new variables: 


K 2 n 2 

n " a 


2 U 

Xk 2 H 2 ' 


( 21 ) 


The coordinates represent the fractional contributions of 
the ordinary energy density, the curvature, the brane ten¬ 
sion, and the dark radiation energy to the total energy 
density, normalized with respect to the Hubble parame¬ 
ter H(t). In addition, we use the variable 7 £ [0,2] in 
order to describe the dynamics of the scalar field <j>. Using 
these variables, the Friedmann constraint (0 becomes: 


f2p + + £lu — 1 . 


( 22 ) 
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Since all of the terms in the sum are non-negative, each 
of the variables fl p , fl k , fl x , flu takes values in the in¬ 
terval [ 0 , 1 ], which means that we obtain a compact state 
space. 

In order to decouple the evolution equation for H © 
from the evolution equations for the density parameters 
fli, we introduce the dimensionless time derivative 

' = IH^d/dt. (23) 

Inserting the Friedmann constraint ( |22| ) in order to keep 
the dimensionality of the state space as low as possible, 
it follows that the dynamics of the open or flat models 
with positive IA are described by the following equations: 


Notice that F+ ^ only occurs for 0 < b 2 < 6 , x/ 3 (&) 

occurs in this sector of state space for b 2 > 2 , and a/ 3 ( 6 ) 
occurs only for b 2 > 4 in this sector of state space. All 
three points move in state space as we vary the parameter 
value b, but independent of the value of b only occur in 
the expanding sector e = + 1 . 

Note that equilibrium points M f °, m®{fl A , Ha) and R° 
are non - hyperbolic; we analyzed their nature using the 
perturbative methods described in Paper 1 ©• 

B. IA > 0 and 3 R > 0 


7 ' = £^ 3^(7 - 2 )[v / 37 + e byj 1 - fl k -fl A -fl x - , Here we wiU use the variables 7 , Q, fl p , H A , “A, “M, 
n' k = e[(3 7 - 2)(1 - fl k ) + 3 7 (H A - n A ) + (4 - 3 7 )H w ]L! fe , where we define 

= e[37(l + fl x — fl A — flk — Qu) + 2 fl k + 4flu)fl A , (24) q =— with D 2 = H 2 + - 3 B C25) 

n' A = e[3'y(n x -n A ~n k -i~n u ) + 2 n k + 4n u }n x , D ’ 

fl' u = e[(3 7 — 4)(1 — flu) + (2 — 3y )fl k + 37(S2a — ^a)]^w - an d the variables with a tilde are the analogues of those 

in (EOl) and (Ell) but normalized with respect to D instead 
The equilibrium points of this dynamical system, their of H The Friedmann equation (g then becomes 
coordinates in state space and their eigenvalues are given — v 


in TABLE I below. Most of the equilibrium points of the 
system ( pd| ) are the same as those obtained in Paper I 
pot : the flat FLRW universe F 2 with stiff matter and 
a(t) = t 1 / 3 ; the Milne universe M 2 with stiff matter and 
a(t) = t; the Milne universe with 7 = 0 (M e °); the flat 
non-general-relativistic model m 2 with 7 = 2 and scale 
factor a(t) = t 1 / 6 , which has been discussed in detail 

in ||]; the flat FLRW universe F+ ^ 3 with 7 = 7 - and 
a{t) = t 2 / b2 ; and a set of universe models X+ 3 (b) with 
7 = | and curvature fl k = 1 — ^ depending on the value 
of b. 

Furthermore we find the de - Sitter model dS 2 with 
7 = 2 and scale factor a(t) = ea;p(y / A/3t), and a set 
of non-general-relativistic critical points with 7 = 0 
which extends both in fl A - and fl\ - direction. This set 
is denoted by m®(fl A ,fl\)-, it contains the flat general 
relativistic FLRW model and the general relativistic de - 
Sitter model with constant energy density p respectively. 
This 2 - parameter set of equilibrium points contains the 
line of non - general - relativistic critical points m°(fl\) 
obtained in our previous paper . 

In addition to these equilibrium points with vanishing 
U , we find the following points with IA > 0: the points 
R 3 and R 2 which have the same metric a(t) = t 1 / 2 as the 
flat the FLRW model with 7 = | (radiation), but have 
7 = 0, 2 respectively and p = k = A -1 = 0. The point 
A+ 3 (b) describes a flat expanding model with 7 = | 
and vanishing brane - tension ( fl\ = 0), but which has 
in general non - vanishing energy density and non-local 
energy density contributions fl p , flu > 0. For 6 = 0, 
this model has maximal positive IA {flu = 1, fl p = 0), 
whereas for 6 2 = 4 the model coincides with the flat 
FLRW model with 7 = 4/3 (fi p = 1, flu = 0). The scale 
factor of this model is proportional to t 1 / 2 for all values 
of b. 


flp + fl A + fl\ + flu — 1 > (26) 

from which it can be seen that the quantities f \ take 
values in [0,1], whereas Q takes values in [—1,1], Hence 
these coordinates define a compact state space. Intro¬ 
ducing the time derivative 

' = D~ 1 d/dt , 

the evolution equation for D decouples from the evolu¬ 
tion equations for flu Inserting the Friedmann constraint 
(^(|), we find that the dynamics of the closed and flat 
models are described by 

7 ' = 1 / 37(7 - 2 ) [\/37 Q + b\Jl-fl\-fl A - fl u ] , 

Q' = [1 — - 7(1 + fl\ — fl A — flu) — 2f2 w ](l — Q 2 ) , 

fl' A = [ 37(1 + fl\ — fl A — flu) + ^flu]Q&A j 
fl’ x = [3 7 (H a - fl A - 1 - flu) + ^fl u ]Qfl x , 

= [(37-4)(1^H w )+3 7 (H a -H a )]QH w . (27) 

The equilibrium points of this system, their coordinates 
in state space and their eigenvalues are given in TABLE 
II. We recover the equilibrium points obtained in the 
previous subsection which corresponded to flat models. 
In particular, for b 2 > 4 we find the equilibrium point 
A / 3 (b), which for b 2 =/= 4 describes a model with non- 
vanishing non-local energy density IA. 

In addition, we find a 2 - parameter set of static {Q = 0) 
models E. These models occur for any fixed values of 
coordinates ( 7 *, 0 , fl A , flu) subject to the constraints 
( p 6 | ) and 

l-^*(l + fl* x -fl* A -fl* u )~2fl*u = 0 (28) 




5 


TABLE I: This table gives the coordinates and eigenvalues of the critical points with U > 0 and 3 R < 0. We have defined 
^ = \/ TP ~ 3 > C = 


Model 

Coordinates 

Eigenvalues 

Fe 

(2, 0,0,0,0) 

(6e + a/6 b, At, 6e, —6e, 2e) 

M° 

(0,1, 0,0,0) 

non-hyperbolic 

Ml 

(2,1,0,0,0) 

2e(3,-2,1,-5,-1) 

ASl 

(2, 0,1, 0,0) 

(6e, —2e, —6e, — 12e, -4e) 

m°(« A ,n A ) forb = 0 a 

(o,o,n A ,n A ,o) 

—2e(3, 1, 0, 0, 2) 

m°(fi A ,ff A ) for b^0 0 

(o,o,n A ,n A ,o) 

(oo 6 , —2e, 0, 0, —4e) 

ml 

( 2 , 0 , 0 , 1 , 0 ) 

2e(3, 5,6, 3,4) 

R° 

( 0 , 0 , 0 , 0 , 1 ) 

non-hyperbolic 

R l 

(2, 0,0,0,1) 

2e(3, 1,2, -4,-1) 

pb 2 /3 c 

X^ 3 (6) d 

(f .0,0, 0,0) 

(f - 3,6 2 - 2,6 2 ,-fo 2 ,6 2 - 4) 

(§,l-£,0,0,0) 

(—1 — ip, —1 + ip, 2, —2, —2) 

A + /3 (fe) e 

(|,0,0,0,1-^) 

(-|(1 + C),2,4, -4,-§(1-0) 


a The eigenvalues of these points can only be evaluated for S2 a + 
7 ^ 1. For n A + = 1 a perturbative analysis has to be carried 

out. 

fe This actually reads lim^^ o(—6e — \/3 byj 1— ^~^ A ) = oo. 

c b 2 e[ 0 , 6 ] 

d b 2 G [2, oo[ 
e b 2 G [4, oo[ 


TABLE II: This table gives the coordinates and eigenvalues of the critical points with U. > 0 and 3 R> 0. We have defined the real 

quantities x = \J 8 ~| &2 ,$ = — ^)\j 3t(1 — — &u) + 2 Clu — 1, a = yj 97 2 (Qa + &u - 1) + 67(1 — 2 Clu) + 4(1 - flu)- 

Notice that a is a real positive quantity within the allowed parameter range. 


Model 

Coordinates 

Eigenvalues 

Fe 

(2, e, 0,0,0) 

(6e + a/ 6b, At, 6e, —6e, 2e) 

dS? 

(2, e, 1,0,0) 

(6e, —2e, —6e, —12e, -At) 

m° (Q.a , ff A ) for b = 0 “ 

(0,e,fj A ,n A ,0) 

—2e(3,1,0, 0,2) 

m?(fi A ,fi A ) forb^O a 

( 0 , e, n A , n A , 0 ) 

(oo 6 , —2e, 0, 0, —4e) 

Re 

(0,e, 0,0,1) 

non-hyperbolic 

ml 

(2, e, 0,1,0) 

2e(3, 5,6, 3, 4) 

R l 

(2, e, 0,0,1) 

2e(3,1,2, -4,-1) 

E 

(7*,o,ni,nx,n&) 

(£,a,0,0, -a) 

p6 2 /3 c 

X^ 3 (6) d 

(t, 1 - 0 , 0 , 0 ) 

( 3 - — 3, b 2 — 2, b 2 , —b 2 ,b 2 — 4) 

(§,-^,0,0,0) 

(^-X,^ + X,-V2b,V2b,V2b) 

A + /3 (fo) e 

(|,l,0,0,l-£) 

(-5(1 + 0,2,4, -4,-i(l-O) 


“The eigenvalues of these points can only be evaluated for £l A + 
Q\ 1. For n A + f2_x = 1 a perturbative analysis has to be carried 
out. 


& This actually reads o(—6e — y/3byj -—) = 00 . 

c b 2 G [0,6] 
d b 2 G [0, 2] 
e b 2 G [4,oo[ 


as well as 

Vr (7* - 2)by/i-n%-ni-C^ = 0 . (29) 

Again, we have already inserted the Friedmann con¬ 
straint (26). Note that the Jacobian of the dynamical 
system (49) is in general not well - defined for 7 = 0 [^lj . 
In that case a perturbative study confirms that the points 


E | 7 *=o are indeed saddle points in state space. 

The models E have the same line - element as the Ein¬ 
stein universe (H = 0 , k = + 1 ), but in general also non¬ 
zero brane tension and positive nonlocal energy density 
U. The quantities a,p,U are constants a*, p* ,U*. 

Requiring positive energy density p *, we find from (|l4|) 
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that the constants a* ,U* must satisfy 


“ 2 V 2 3 ' 


(30) 


Notice that the constraint (H|) conies from the j' - 
equation in (p7j), hence from the dynamics of the scalar 
field, and does not occur in the studies of a similar model 
with constant equation of state, i.e. without a dynamical 
equation of state parameter. 

From (H|) we can see that the static models E can oc¬ 
cur for any values of 7 £ [0, 2]. Equation (|2f]) is however 
in general only satisfied for all values of 7 for the vacuum 
models (H p = 0) [ p2| . 

Note that in contrast to dynamical models, it is not 
possible to have fi p = 0 and fl\ ^ 0 when considering 
static models. This is due to the fact that for static 
models a = a*, p = p* are constants and therefore H p = 
0 implies p* = 0, which leads to Ha = 0. This has the 
interesting consequence that there exist no vacuum static 
solutions in the A = 0-subset. Only when allowing for 
the additional degree of freedom Ha do we find that ( |28| ) 
can be satisfied for H p = 0. 

We emphasize that the static model denoted by E 1 / 3 in 
our previous paper [|lo]] is in that sense unphysical, since 
that model found in the A = 0 case has Q = 0, H p = 0 
and fl\ = 1 which can only be satisfied in this extended 
scenario allowing for non-zero cosmological constant A. 

Hence the only static vacuum model occurring in this 
sector is the model with Ha = 0, which means that 
flu = 1/2, Ha = 1/2. This is a line of equilibrium points 
extending in 7 - direction. The other static models are 
the ones with 7 = 0 or 7 = 2. The former are character¬ 
ized by flu = 1/2, the latter have flu = 2 + 3(Ha — Ha). 

Note that the only static models with Ha = 0 are the 
ones with 7 = 0, flu = 1/2 , Ha + H p = 1/2. 

For constant potential (b = 0) equation ( p9| ) is auto¬ 
matically fulfilled for all values of 7*, fl* K , H^, H^. Thus 
for b = 0 E degenerates into a 3 - dimensional surface 
containing static models for all values of 7 £ [0,2]. 

Notice that in the Ha = 0-subset much stronger con¬ 
straints must again be satisfied; it can be shown that here 
7 £ [0, |] is required even for b = 0. 

Although these sets of equilibrium points form geomet¬ 
rically interesting objects in state space, they are not of 
interest for our stability analysis, since all these equilib¬ 
rium points are unstable saddle points. 


and the variables Hj are defined as in in (^) and ( |2l| ) 
but normalized with respect to N instead of H. Using 
these variables, the Friedmann constraint reads 

H p + flk + Ha + Ha = 1 (32) 

As before, all the terms in that sum are non - negative by 
definition, hence the variables H p , Hfc, Ha, Ha take values 
in the interval [0,1]. Since 7 £ [0,2] and Z £ [—1,1], we 
find that the state space is again compact. 

Furthermore, introducing the time derivative 

' = N~ 1 d/dt , (33) 

the evolution of N decouples from the rest of the vari¬ 
ables. We obtain 

i = - 2) [\/37 Z + b\Jl - fl x - fl k - H a ] , 

Z’ = —[^7(1 + H A — H a ) + (1 — ^ 7 )Hfe — 2](1 — Z 2 ) , 

fl'k = [(3 7 — 2)(1 — flk) + 3 7 (H a — H A )]^Hfc , 

Ha = [3 7 (1 - H fc - H a + H a ) + 2H fc ]ZH A , 

H'a = [2H fc - 3 7 (1 + H fc + H A - fl\)]Zfl x . (34) 

The equilibrium points of this system, their coordinates 
in state space and their eigenvalues are given in TABLE 
III. Again we recover the models with U = 0 that we 
have obtained in the previous two sections. Notice that 
the point A+ 3 (b) now only occurs for b 2 £ [0,4]. This 
is due to the fact that the point represents a model with 
nonlocal energy density depending on the value of the 
parameter b. For b 2 < 4, the point corresponds to a 
model with negative nonlocal energy density U, whereas 
for b 2 > 4 it represents a model with positive U. For 
b 2 = 4, the point coincides with the expanding FLRW 
model F+ ,/-3 | b 2 =4 with U = 0 and 7 = |. Thus the 

point At{_ 3 (b) moving in state space is leaving the sector 
describing models with U < 0 and entering the sector of 
models with U > 0 when the parameter value b 2 = 4 (6 = 
-2). 

In addition to these points, we find the set S of static 
models (H = 0) whose coordinates (7*, 0, Hj£, fl\, H]() 
satisfy the constraints 

2 - | 7 *(1 + Ha - H^ — fiJO — HJ = 0 (35) 


and 


C. U < 0 and Z R < 0 

In this case we obtain a compact state space by intro¬ 
ducing the dynamical variables Z. H p , Hfc, Ha, Ha, where 


Vr (7 * - 2)5^/i-H* -fl\~fl\ = 0, (36) 

where we have already inserted the Friedmann constraint 

Requiring non - negative energy density p *, we also find 
from (|l4|) that the constants a*,U* must satisfy 



N 2 = H 2 — 


(31) 


-w > E-A 

- 2 v 3 



2 U 
A/? ’ 


(37) 
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TABLE III: This table gives the coordinates and eigenvalues of the critical points with U < 0 and 3 R < 0. We have defined ip = 
^-3, 5 = V64 — 15b 2 ,7] = v/26(| - i)^^ _ fi A _ fi fc ) + _ 2 , T = y/9^(Q. A + fi fc - 1) + 6 7 (2 - fi fc ) + 4(1 - tl k ). 

Notice that T is a real positive quantity within the allowed range of variables. 


Model 

Coordinates 

Eigenvalues 

F 2 

(2, e, 0,0,0) 

( 6 e + 1 / 66 , 2e, 4e, 6e, —6e) 

M° 

(0,6,1,0,0) 

non-hyperbolic 

M 2 

(2, e, 1,0,0) 

2e(3,-1,-2,1,-5) 

dS 2 

( 2 ,6, 0 , 1 , 0 ) 

(6e, —4e, —2e, -6e,-12e) 

m°(n A , ri\) for b — 0 " 

( 0 , e, 0 , n A , ri\) 

—2e(3, 2,1,0, 0) 

m°{Q. A ,rix) forb^ 0 “ 

(0,e,0,n A ,fiA) 

(oo\—4e, -2e,0,0) 

m? 

( 2 , e, 0 , 0 , 1 ) 

2e(3,4, 5,6, 3) 

S 

( 7 *,o 

(?7,r,o,-r,o) 

F+ /3 

(t* 1 . 0 , 0 , 0 ) 

(4y — 3, fo 2 — 4, b 2 — 2, b 2 , —b 2 ) 

X 3/3 (fc) c 

(§,l,l-£,0,0) 

(—1 — ip, —2, —1 + ip, 2, —2) 

A d/3 (6) d 

( 3 ,- 5 , 0 , 0 , 0 ) 

{\{b-£,),\(b + ii),-b,-2b,2b) 


“The eigenvalues of these points can only be evaluated for fl,\ + 
Clx 7 = 1 . For T fix = 1 a perturbative analysis has to be carried 
out. 

^This actually reads o(—6e — \/3byj 1 ~^ A ~ n v .) _ 00 _ 

c b 2 S [2, oo[ 
d fe 2 e [o, 4] 


It is important to note that these static models in general 
have negative curvature (k = —1). 

We can see from ([f|) and ([l2|) that there are no static 
vacuum models. The only static models for i) / 0 are 
the ones with 7* = 2 and | + |(fi A ~ ^a)- These 
models occur both in the Cl A = 0-subset and outside of 
that region. Again this is a 2-parameter surface. 

For 6 = 0, this 2 - dimensional set degenerates into a 
3 - dimensional surface; in this special case we find that 
static models occur for all values of 7* > 2/3. 

Notice that the constraints on these models, in contrast 
to the constraints on the models E found in the previous 
section, do not change when we restrict ourselves to the 
A = 0-subset: for b = 0 there are static models for all 
7 > 2/3 even for fi A = 0. 

Dynamically these models are again not very interest¬ 
ing, since they all represent saddle points in state space. 


D. U < 0 and 3 R > 0 


As explained in Q, we need to take into account that 
the Friedmann equation © now has two non - positive 
terms. It turns out that we can obtain a compact state 
space introducing the dimensionless dynamical variables 
W, n p , fl a, ri\, flu, where 



P 2 = H 2 


V 2 U 
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and the variables with a hat are defined as those in ( |2(i| ) 
and ( pl| ) but normalized with respect to P instead of H. 
The Friedmann constraint now reduces to 


rip + ri A -f ri\ — 1. 


It is important to note that f In does not appear in the 
Friedmann equation ( |38| ) • Indeed, it can be seen from its 
definition that f lu is negative and belongs to the interval 
[—1,0]. As before, we can see from ( |38|) that Cl p , Cl A , (l\ £ 
[0,1]. Together with W £ [—1,1] and 7 £ [0,2], we find 
that these variables define a compact state space. 

Using the time derivative 

' = P~ 1 d/dt, (39) 

we obtain 

i = - 2 )[\/37 W + bsjl - ri x - U A ] , 

W' = [1 — ^7(1 + Ua ^ U A )](1 — W 2 ) — flu , 
fi' A = 37(1 + fix - ri A )wri A , 
ri'x = 3 7 (u A -riA- i)wrix, 

ri' u = [3 7 (1 + Ua-U a )- 4]IUU w . (40) 

The equilibrium points of this system, their coordinates 
in state space and their eigenvalues are given in TABLE 
IV. Again, we recover the equilibrium points that repre¬ 
sent flat models obtained in the previous subsection. In 
addition, we find another set of equilibrium points rep¬ 
resenting static models denoted by E. 

These models correspond to the Einstein-universe - 
like models found in subsection B, since they also have 
the same line element as the Einstein universe (H = 
0, k = +1) and have non-zero brane tension and non¬ 
local energy density U, except now U < 0. 

Their coordinates (7*, 0, (l* A , satisfy the con¬ 

straints 

1 - | 7 *(1 + ^a - ^a) - = 0 


(38) 


(41) 














TABLE IV: This table gives the coordinates and eigenvalues of the critical points with U < 0 and 3 R > 0. We have defined 

X = \J an d £ = \/6?-M56^, 5 = \/6b\/ 1 — Da — Da, fi = \J 97 2 Da + 4. Notice that /3 is a real positive quantity within 
the allowed range of variables. 


Model 

Coordinates 

Eigenvalues 

F i 

(2, e, 0,0,0) 

(6e + a/ 6b, 4e, 6e, — 6e, 2e) 

dS 2 

(2, e, 1,0,0) 

(6e, — 2e, -6e, -12e, -4e) 

m/(f2A, Da) for b = 0 a 

( 0 , e, D a , Da, 0 ) 

-2e(3,l,0,0,2) 

m°(D A ,D A ) forb^0 a 

(0, e, D a , Da, 0) 

(oo 6 , —2e, 0, 0, —4e) 

m? 

(2, e, 0,1,0) 

2e(3, 5, 6, 3,4) 

E 

( 7 *, 0, D^, Dy, Dy) 

(S,P, 0,-/3,0) 

F + /3 

(t. 1 . 0 . o,°) 

(^- — 3, 6 2 — 2, 6 2 , —b 2 ,b 2 — 4) 

X+ /3 (fo) c 

(§,-75,0,0,0) 

(72 - X, 75 + X, -V2b, V2b, V2b ) 

A+ /3 (&) d 

( 5 ,-1,0,0,$-i) 

(z(b-0,z(b + 0,~b,-2b,2b) 

“The eigenvalues of these points can only be evaluated for + 


^ 1. For Oa + ^a — 1 a perturbative analysis has to be carried 


out. 



M il is actually reads tm,_ 0 (— 6s 

c b 2 e [ 0 , 2 ] 
d b 2 e [ 0 , 4] 

V3b^J 1 - W - W ) = qo . 


and 

that in the A 

= 0 -subset the general-relativistic models 


Vr ( 7 * - = 0 . 


(42) 


where we have again inserted the Friedmann constraint 

We can see from © that there exist no vacuum solu¬ 
tions and no solutions for 7 * = 0. Hence for a non-flat 
potential the only static models are the ones with 7 * = 2 
and fly = — 2 — 3 ( 0 ^ — ^a)- particular there are no 
static models for 0 / = 0 . 

For b = 0, we find static models for all 7 > 1/3; the 
ones in the Q/ = 0-subset occur for 7 £ [1/3,4/3]. 


E. Qualitative Analysis 

We use the results obtained in the previous subsections 
in order to determine the dynamical character of the equi¬ 
librium points found above. We summarize the results in 
TABLE V below. For all values of b the past attrac¬ 
tors of the FLRW models are the BDL models m\ with 
7 = 2 ( V « (j) 2 ) and the set of non-general relativistic 
models (Qa,Qa) with 7 = 0 (</ 2 << V). For a non¬ 
flat potential (6 ^ 0) the contracting BDL model m 2 __ 
with 7 = 2 is the unique future attractor. In particular 
the expanding deSitter model is not a future attractor in 
the presence of a scalar field with exponential potential. 
This stands in strong contrast to the results obtained in 
|| for the scalar field - free case where the expanding de - 
Sitter model dS + is a future attractor for all values of 7. 
In the special case of a flat potential (6 = 0) we find that 
the models m[]_(Q a,^a) including the de-Sitter model 
with 7 = 0 form another set of future attractors. Note 


F+ / 3 and X+ 3 (b) are also future attractors. To be pre- 

^2 j o 

cise, the flat FLRW model F + ' is an attractor for 6 2 < 2 
and the model X+ 3 (b) for b 2 > 2. 

Notice that in General Relativity, the scenario that 
we are describing here (matter described by a dynami¬ 
cal scalar field with exponential potential) only admits a 
static universe in the special case of a constant potential 
(b = 0). This model had the same line element as the 
Einstein universe {H = 0, k = +1) and only occurred at 
7 = 2/3. 

If we allow for non - zero brane tension but neglect the 
bulk effects (U = 0, see Paper I ( fid ) )), this condition is 
relaxed: we find static universe models characterized by 
FI = 0, k = +1 for all values of 7 e] 1/3, 2/3] for b = 0. 
For a non - flat potential (b ^ 0) there are still no physical 
static equilibrium points ]^3| . 

In this paper taking into account the bulk effects as 
well as allowing for a non - zero cosmological constant A, 
we find that there are not only static Einstein universe 
like models (k = + 1 ), but also static saddle points that 
are flat {k = 0 ) or even negatively curved (k = — 1 ). 

For b = 0 there are static models for all values of 7 
even in the A = 0-subset, which shows that these are 
purely due to the bulk effects. Allowing for A / 0 we 
find that the constraints on the static models are further 
relaxed. We then find that Einstein static models occur 
for all values of 7 when U > 0 and for all 7 > 1 /3 when 
U < 0. Flat and open static models occur for all 7 > 2/3 
but only if U <0. 

More interestingly we also find static models for b 7 ^ 0. 
For A = 0 we find Einstein static models with 7 = 0 in 
the U > 0-sector and for 7 = 2 there are open and flat 
models in the U < 0-sector. 
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TABLE V: Dynamical character of the critical points in the FLRW case for A = 0. 


Model 

b = 0 

0 < b 2 < 2 

b 2 =2 

2 < b 2 < 4 

b 2 = 4 

4 < b 2 < 6 

b 2 — 6 

b 2 > 6 


sink 

saddle 

saddle 

saddle 

saddle 

saddle 

saddle 

saddle 

m° (Q x ) 

source 

source 

source 

source 

source 

source 

source 

source 

E 

saddle 

saddle 

saddle 

saddle 

saddle 

saddle 

saddle 

saddle 

S 

saddle 

saddle 

saddle 

saddle 

saddle 

saddle 

saddle 

saddle 

F? 

saddle 

saddle 

saddle 

saddle 

saddle 

saddle 

saddle 

saddle 

mi 

source 

source 

source 

source 

source 

source 

source 

source 

m 2 _ 

sink 

sink 

sink 

sink 

sink 

sink 

sink 

sink 

R° 

saddle 

saddle 

saddle 

saddle 

saddle 

saddle 

saddle 

saddle 

R? 

saddle 

saddle 

saddle 

saddle 

saddle 

saddle 

saddle 

saddle 

Fi 2/3 

sink 

sink 

saddle “ 

saddle 

saddle 

saddle 

saddle 

- 

CO 

saddle 

saddle 

saddle “ 

sink 

sink 

sink 

sink 

sink 

A l /3 (b) b 

saddle 

saddle 

saddle 0 

saddle 

saddle 

saddle 

saddle 

saddle 

“Notice that this point 

is an attractor for all general relativistic 





open or flat 

models, while it is a repeller for all closed general 





relativistic models. 








6 Notice that this model has negative non- 

-local energy density IA 





for b 2 <4, U 

= 0 for b 2 = 

= 4 and U > 0 for b 2 > 4 








TABLE VI: Dynamical character of the critical points 

in the FLRW case for general A. 


Model 

O 

II 

-O 

0 < b 2 < 2 

b 2 = 2 

2 < b 2 < 4 

b 2 = 4 

4 < b 2 < 6 

b 2 = 6 

b 2 > 6 

hi^HajHa) 

sink 

saddle 

saddle 

saddle 

saddle 

saddle 

saddle 

saddle 

m° (fi A , H A ) 

source 

source 

source 

source 

source 

source 

source 

source 

E 

saddle 

saddle 

saddle 

saddle 

saddle 

saddle 

saddle 

saddle 

S 

saddle 

saddle 

saddle 

saddle 

saddle 

saddle 

saddle 

saddle 

F 2 

saddle 

saddle 

saddle 

saddle 

saddle 

saddle 

saddle 

saddle 

dS 2 

saddle 

saddle 

saddle 

saddle 

saddle 

saddle 

saddle 

saddle 

m+ 

source 

source 

source 

source 

source 

source 

source 

source 

m 2 _ 

sink 

sink 

sink 

sink 

sink 

sink 

sink 

sink 

R° 

saddle 

saddle 

saddle 

saddle 

saddle 

saddle 

saddle 

saddle 

R? 

saddle 

saddle 

saddle 

saddle 

saddle 

saddle 

saddle 

saddle 

CO 

^ + 

sink 

saddle 

saddle 

saddle 

saddle 

saddle 

saddle 

- 

X 2/3 (6) 

saddle 

saddle 

saddle 

saddle 

saddle 

saddle 

saddle 

saddle 

A+ /3 (&)“ 

saddle 

saddle 

saddle 

saddle 

saddle 

saddle 

saddle 

saddle 


“Notice that this model has negative non-local energy density U 
for b 2 < 4, U = 0 for b 2 = 4 and U > 0 for b 2 > 4 


When allowing for a cosmological constant, we find 
that the Einstein static universe occurs for all values of 
7 when U > 0 and for 7 = 2 when U < 0. The open 
and flat models still only occur for 7 = 2 and only in the 
U < 0-sector. 


Altogether, the conditions for allowing for static mod¬ 
els in General Relativity are changed dramatically when 
considering the terms corresponding to brane tension and 
non - local energy density, and are further relaxed when 
including a cosmological constant. Instead of finding only 
the static Einstein universe with 7 = |, and this only in 
the special case of a flat potential, we now find static 
models for all values of b and all 7 € [0, 2]. 


IV. BIANCHI I MODELS WITH AN 
EXPONENTIAL POTENTIALS 

We now turn our attention to the dynamics of Bianchi 
type I models in the brane - world scenario with exponen¬ 
tial potential. This class of models is characterized by a 
metric of the form 

3 

ds 2 = — dt 2 + ^ A 2 a {t){dx a )‘ 2 . (43) 

Q! = l 

Again, we follow very closely the analysis done in f|. the 
only difference being that we have the additional equa¬ 
tion ( 0 ) describing the dynamics of the scalar field. 

As we have discussed above, the non - zero contribu¬ 
tions from the five - dimensional Weyl tensor are U and 
V a b but since the second one has no evolution equation 
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we will assume for simplicity that V a bO' ab = 0 . 

When introducing appropriate variables, we have to 
consider two different cases: (A) U > 0; (B) U < 0. 


A. U > 0 


This case is formally very similar to subsection III A. 


We obtain a compact state space using the dimensionless 
variables 7 , f 2 p , 11 a, II A , , where 


= 


3 h 2 ' 


(44) 


defines a normalized variable for the shear contribution, 
and the remaining variables are defined by and ( 21 ). 
The Friedmann constraint 0 ) reduces to 


With respect to the time derivative ' = N 1 d/dt , the 
dynamical system becomes 

7' = - 2 ) W&yZ T b\J\ &\ fl A ^<t] 1 

Z' = -[(|^ - 2) -+- (3 - + |^(^A - 0 A )](1 - ^ 2 ) , 

n' A = [(6 - 37)n ff + 37(1 + n A - ci A )]zn A , 

= [(37-6)(i-n <T ) + 37(n A -n A )]zn <T , 

= [(6-3 7 )n CT + 37(0 A -f) A -l)]^ A . (49) 

The equilibrium points of this system, their coordinates 
in state space and their eigenvalues are given in TABLE 
VII. In addition to the anisotropic Kasner models we 
also obtain the set of saddle points nE whose coordinates 
( 7 *,0,fiX,f^,O*) satisfy 

2 - §7*(1 + - O* - HI) - 30* = 0 (50) 


£lp + + Qu — 1 • (45) and 


Using the time derivative 


Vrh* - 2)byJi-cr A -n* x -n %]. ( 51 ) 


' = \H\- x d/dt 

the system of dynamical equations is given by 


7 ' = £1/37(7 - 2 ) [a /37 + e by/l ~ - n A - - Slu\ 

fl' A = e[3 7 (l + f2 A — fl A ) + (6 — 3 7 )f2 CT + (4 — 3 7 )f2^]fI A 
K = £ [(3 7 — 6)(1 — 12^) + 3 7 (O a — fI A ) + (4 — 3 7 )f2^]f2, 
= e[37(^A - - 1) + (6 - 37 )^ 0 . + (4 - 37)fI M ]fi A 

% = e [(3 7 — 4)(1 — flu) + (6 — 3 7 )fl, 


Furthermore the constants U*, a* 2 must satisfy the con¬ 
dition 


-U* > 


A k 2 


(c 


•A) 


(52) 


The equilibrium points of this system, their coordinates 
in state space and their eigenvalues are given in TABLE 
VI. The only new equilibrium points that we obtain in 
this subsection are the anisotropic models K. These have 
been identified in || as the general relativistic vacuum 
Kasner models with line element (O) where 


in order for the energy density p* to be non-negative. 

We can see from (50) and (4^) that there are vacuum 
’solutions with S2* = 2/3 for al values of 7 *. Furthermore 
we find non-vacuum models for 7 * = 0 and 7 * = 2 . 
37 (S 2 A — Ha)]^w The set of static models nE forms a 2-parameter set of 
( 46 'equilibrium points which for 5 = 0 degenerates into a 3- 
dimensional surface containing static equilibrium points 
for all values of 7 *. 

Note that the Cl A = 0 - subset only contains static mod¬ 
els for 7 * G [ 0 , |] if 6 yA 0 and for 7 * e [ 0 , |] if b = 0 . 

For each fixed value of 7 * the models nE are the same 
as the ones found in [Bl with scale factors 


A a {t) = e 1 


q a t 


3 3 

A a = t 2pa and 0 Pa = 0 pI = f ■ (47) 


a— 1 


a—1 


B. U < 0 


In analogy to subsection IIIC , we consider the vari¬ 
ables 7 , Z , fL, VL\, ff) A , where Cl a is defined as f l a in 
equation (|4J), but normalized with respect to N defined 
in and the remaining variables are as in subsec¬ 
tion HI C. With these variables the Friedmann constraint 


(fh|) reads 


flp + + f i a + fl A — 1 . (48) 


where q a are constants which satisfy 
3 3 

0 = 0 and = 2cr * 2 ■ 

a=l a=l 

We repeat that although H = 0, the scale factors A a (t) 
are in general not static. It is only the overall scale factor 
J )/ =1 A a (t) = 1 which remains constant. The matter can 
still expand or contract along the principal shear axis, 
hence forming a pancake singularity if one q a is negative 
and a cigar type singularity if two of them are negative. 

Note that the Jacobian of the dynamical system (|49| ) 
is in general not well defined for 7 * = 0 (2l|] . In that case 
a perturbative study confirms that the points nE | 7 *=o 
are also saddle points in state space. 












11 


TABLE VII: This table gives the coordinates and eigenvalues of the critical points with U > 0. We have defined £ = , /1§ — 15. 


Model 

Coordinates 

Eigenvalues 

Ke 

(0,0,1,0,0) 

non-hyperbolic 

k ?(n CT ) 

(2,o,n<r,o,o) 

(6e + V66\/l — H<r, 6e, 0, —6e, 2e) 

dS? 

(2,1,0, 0,0) 

(6e, —6e, —6e, — 12e, —4e) 

m°(fl A , fi A ) for 6 = 0" 

(0, 0 A , 0, 0 A , 0) 

—2e(3, 0,3, 0,2) 

m°(fi A , fi A ) for 6^0° 

(o,n A ,o,n A ,o) 

(oo 6 , 0, —6e, 0, —4e) 

m? 

(2, 0,0,1,0) 

2e(3,6,3,3,4) 

R° 

(0,0, 0,0,1) 

non-hyperbolic 

R? 

(2, 0,0, 0,1) 

2e(3,2,-1,-4,-1) 

pb 2 /3 c 

Af" 3 (6) d 

(t. 0 . 0 ,0,0) 

(£ - 3,6 2 ,6 2 - 6,-6 2 ,6 2 - 4) 

(S,0,0,0,l-£) 

(-|(1 + 0,4, -2, -4,-|(l-0) 

a The eigenvalues of these points c 

:an only be evaluated for f^A + 


^ 1. For Qa + ^a = 1 a perturbative analysis has to be carried 


out. 



6 This actually reads o(—6e 

c b 2 e [ 0 , 6] 
d b 2 e [4,oo[ 

73 = 00 . 


TABLE VIII: This table gives 

the coordinates and eigenvalues of the critical points with U < 0. We have defined £ = 

V64 - 156 2 , 0 = V2b(f - ±)\/37(l - 0 A - fi*) +30,-2, ip = 79 7 2 (fl A + 
real positive quantity within the allowed range of variables. 

— 1) — 187fl CT + 12y + 4. Notice that ip is a 

Model 

Coordinates 

Eigenvalues 

Ke 

(0,e, 0,1,0) 

non-hyperbolic 

K?(0,) 

(2, e, 0,ft«,,0) 

(6e + 766v 7 ! - 2e, 6e, 0, -6e) 

dS? 

(2, e, 1,0,0) 

(6e, — 4e, — 6e, — 6e, — 12e) 

m?(n A ,n A ) forb = 0‘ 

(0,e,ft Al 0,ft A ) 

— 2e(3, 2, 0, 3, 0) 

m°(O\,O0 f or b ¥=0 a 

(o, e ,n A ,o,n A ) 

(oo 6 , — 4e, 0 , — 6e, 0 ) 

m? 

( 2 , + 0 , 0 , 1 ) 

2e(3,4,6,3,3) 

nE 

( 7 *,o,nx,n;,nj) 

{(/),<p,0,-cp,0) 

p6 2 /3 c 

A^ 3 (6) d 

(T, 1 , O’ 0 ’ 0 ) 

(£ - 3,6 2 - 4,6 2 ,6 2 - 6,-b 2 ) 

(5,-5, 0,0,0) 

(l(&- C),i (6 + 0,-26, 6, 26) 

“The eigenvalues of these points c 

;an only be evaluated for ^a + 


7 ^ 1. For = 1 a perturbative analysis has to be carried 


out. 



6 This actually reads o(—6e 

c b 2 e [ 0 , 6] 
d b 2 e [0, 4] 

Vsb^ 1_nA ~^ x ) ~ 00 . 



C. Qualitative Analysis 


All these models are isotropic. 


We summarize the dynamical character of the equi¬ 
librium points obtained in this section in TABLE VIII 
below. The main question we want to address here is 
whether the initial singularity is isotropic or not when 
restricting our analysis to Bianchi I models. We obtain 
the following result. 

For 6 = 0 we find that the only equilibrium points 
of the dynamical system (46, fhl) that are sources and 


therefore present stable initial configurations are the ex¬ 
panding non-general relativistic BDL model with 7 = 2 
(m+) and the line of the flat and collapsing models with 
7 = 0 (to^_(Qa, fi(A))) including the flat FLRW model 
and the maximally non-general relativistic BDL model. 


If b < 0, the dynamical system possesses in addition 
to the sources above one further source denoted by . 
This equilibrium point represents the anisotropic Kasner 
model with 7 = 0 . 


This means that if b = 0, the initial singularity must 
be isotropic, since the anisotropic models are not stable 
at early times. If on the other hand the potential V((f >) 
is not flat, i.e. 6 7 0 , we find that the initial singularity 
could be anisotropic, since the anisotropic Kasner model 
is stable at early times. 























12 


TABLE IX: Dynamical character of the critical points in the Bianchi I case for A = 0. 


Model 

6 = 0 

0 < 6 2 < 4 

6 2 = 4 

4 < 6 2 < 6 

6 2 = 6 

6 2 > 6 

K° + 

saddle 

source 

source 

source 

source 

source 

K° 

saddle 

saddle 

saddle 

saddle 

saddle 

saddle 

K 2 (f2 CT ) 

saddle 

saddle 

saddle 

saddle 

saddle 

saddle 

m^fU) 

sink 

saddle 

saddle 

saddle 

saddle 

saddle 

m° (fix) 

source 

source 

source 

source 

source 

source 

m + 

source 

source 

source 

source 

source 

source 

m 2 _ 

sink 

sink 

sink 

sink 

sink 

sink 

Re 

saddle 

saddle 

saddle 

saddle 

saddle 

saddle 

Re 

saddle 

saddle 

saddle 

saddle 

saddle 

saddle 

nE 

saddle 

saddle 

saddle 

saddle 

saddle 

saddle 

F+ /3 

sink 

sink 

saddle “ 

saddle 

saddle 

- 

A l /3 (b) b 

saddle 

saddle 

saddle “ 

sink 

sink 

sink 

“Notice that this point is an attractor for all general relativistic 
open or flat models, while it is a repeller for all general relativistic 
closed models. 

6 Notice that this model has negative non-local energy density IA 
for b 2 < 4, U = 0 for b 2 = 4 and U > 0 for b 2 > 4 






TABLE X: Dynamical character of the critical points 

in the Bianchi I case 

for general A. 


Model 

6 = 0 

0 < 6 2 < 4 

6 2 = 4 

4 < 6 2 < 6 

6 2 =6 

6 2 > 6 

K| 

saddle 

source 

source 

source 

source 

source 

K°_ 

saddle 

saddle 

saddle 

saddle 

saddle 

saddle 

K?(fi„) 

saddle 

saddle 

saddle 

saddle 

saddle 

saddle 

dS 2 

saddle 

saddle 

saddle 

saddle 

saddle 

saddle 

m+(n A ,nA) 

sink 

saddle 

saddle 

saddle 

saddle 

saddle 

m° (Da, VL\) 

source 

source 

source 

source 

source 

source 

m + 

source 

source 

source 

source 

source 

source 

m 2 

sink 

sink 

sink 

sink 

sink 

sink 

R° 

saddle 

saddle 

saddle 

saddle 

saddle 

saddle 

R? 

saddle 

saddle 

saddle 

saddle 

saddle 

saddle 

nE 

saddle 

saddle 

saddle 

saddle 

saddle 

saddle 

F+ 73 

sink 

saddle 

saddle 

saddle 

saddle 

- 

A+ /3 (b)“ 

saddle 

saddle 

saddle 

saddle 

saddle 

saddle 


“Notice that this model has negative non-local energy density U 
for b 2 < 4, U = 0 for b 2 = 4 and U > 0 for b 2 > 4 


V. DISCUSSION AND CONCLUSION 

We end this paper with a comparison of our results 
with work previously done in this area. In particular, we 
focus on two issues: Firstly we discuss the occurrence and 
stability of equilibrium points in the dynamical systems 
analysis of FLRW and Bianchi I models with and without 
an exponential potential and compare our work to that 
done by Campos and Sopuerta Q. Secondly we discuss 
the issue of isotropization in the past; here we mainly 
comment on the results recently obtained by Coley et al 
in (|| H|). 

Addressing the first issue let us state the differences 
between this work and || jGj where the scalar field free 
analog of the scenario was studied. Formally, the new 
features of our work is contained in an additional evo¬ 
lution equation for 7 © thus enlarging the dynamical 


system analyzed with respect to the one studied in [^). 
This means in particular that we should expect to obtain 
no more than the equilibrium points already obtained in 
H. The stability of these equilibrium points on the other 
hand is expected to be altered since the additional dy¬ 
namical equation yields an additional eigenvalue in the 
dynamical systems analysis which can potentially desta¬ 
bilize the equilibrium points. In fact we found that we in 
general only recover the equilibrium points that Campos 
et al found for any linear barotropic equation of state 
for the special values of 7 = 0,2 which correspond to 
(j) « V and V « respectively. 

It is worthwhile mentioning that the ’’new” equilib¬ 
rium points F+ 73 , X+ 3 (b) and A+ 3 (b) which are mov¬ 
ing in state space as the steepness of the potential (char¬ 
acterized by the value of b) is increased are in fact not 
new. The first point simply corresponds to the flat FLRW 
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model with equation of state parameter 7 depending on 
the value of 6 . More interestingly, the last two corre¬ 
spond to the bifurcations occurring for 7 = |, | in [|j : 

for each value of b the model xV 3 (b) corresponds to a 
point of the line of general-relativistic, non-static equi¬ 
librium points which occurs in jb| for 7 = § (see figure 

4 in that paper), and similarly the model A+ 3 (b) corre¬ 
sponds for each value of b to one of the points on the 
line of equilibrium points joining the models denoted by 
F + and R + for 7 = | in || (see figure 2 in that paper). 
There the stability of the equilibrium points forming the 
bifurcations was not discussed in much detail. We find 
that for sufficiently steep potential ( 6 2 > 4) the model 
A+ 3 (b) is a stable sink in the A = 0-subset of the state 
space of the Bianclri I models and the A = k = 0-subset 
of the state space of the FLRW models. This means that 
in that case this point which has scale factor a(t ) = f 1//2 
and non-vanishing positive non-local energy density IA is 
a future attractor with the equation of state of pure ra¬ 
diation. It is the only non - collapsing future attractor, 
but it is unstable for A, k 7 ^ 0. 

The second important issue we want to address here 
is the issue of isotropization. In Jl8| it was claimed that 
the initial singularity in the brane-world scenario is nec¬ 
essarily isotropic. It was shown that for 7 > 1 the non¬ 
general relativistic BDL-model m+ j25| with scale factor 
a(t) = t 1 / 37 is a source in the state space of all Bianchi 
IX models, and it was also claimed that this is the only 
source with physically relevant values of 7(7 > 1). In 
this paper we have considered the more general situation 
of 7 € [0, 2], We have found all equilibrium points of the 
state space and identified the past attractors without im¬ 
posing any constraints on 7 . We confirm that the BDL 


model with 7 = 2 is a generic past attractor in both the 
FLRW case and the Bianchi I case. We want to point 
out however, that in the analysis of the Bianchi I models 
the anisotropic Kasner model with 7 = 0 (if?) is another 
past attractor unless the potential is flat (6 = 0). This 
model has the equation of state p = —p and corresponds 
to a slow rolling scalar field (</><< V). In particular no 
energy conditions are violated |2(3| l . In summary, we find 
that if we adopt the assumption on the equation of state 
in the early universe 7 > 1 , the initial singularity must 
be isotropic. The benefit of our approach is that we first 
find all configurations that are stable in the past in a 
transparent and complete analysis. We can then exclude 
certain ones on physical grounds. 

We conclude with the following remark on the future 
attractors in this scenario. In [l0| we explained that the 

general relativistic models F+ and xV 3 (b) correspond 
to the equilibrium points obtained in ]Tl|| and |l^| in the 
general relativistic analysis, and found that the stability 
of these models is not altered in the brane-world exten¬ 
sion with vanishing non-local energy density IA. In this 
paper we find that these models are also stable against 
perturbations in IA. It is the perturbations in A that grow 
and hence destabilize these models. This confirms the in¬ 
tuitive idea that the cosmological constant dominates at 
late times since we showed that allowing for non-zero A 
destabilizes all expanding future attractors with = 0 . 
It does not destabilize the contracting non - general rel¬ 
ativistic future attractor m 2 , since at high energies the 
A term is negligible with respect to the p 2 term corre¬ 
sponding to f 1\. Hence the re-collapsing BDL model 
with stiff matter (7 = 2 ) is the unique future attractor 
for the models discussed in this paper. 
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